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We investigate the positivity of the Euclidean path integral measure for low-energy 
modes in dense fermionic matter. We show that the sign problem usually associated 
with fermions is absent if one considers only low-energy degrees of freedom. We 
describe a method for simulating dense QCD on the lattice and give a proof using 
rigorous inequalities that the color-flavor locked (CFL) phase is the true vacuum 
of three flavor, massless QCD. 



1. Introduction 

Euclidean quantum chromodynamics (QCD) with a non-zero chemical po- 
tential has a complex measure, which has made lattice simulation partic- 
ularly difficult^. (Lattice simulation of the QCD phase boundary at finite 
density^''^ has been a topic of recent interest.) This problem is often re- 
ferred to as the sign problem, because, by appropriately grouping terms, 
quantities such as the partition function can be written as a sum over real, 
but potentially negative, terms. (That this grouping can be accomplished 
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is in many systems a consequence of a discrete symmetry such as parity or 
time-reversal invariance.) Indeterminate signs are enough to preclude use 
of importance sampling, the main technique for speeding up Monte Carlo 
integration. It is important to note that while the sign problem often arises 
in systems of fermions, it is neither inevitable nor inescapable. For exam- 
ple, in QCD at zero chemical potential and in the Hubbard model at half 
filling one can organize the sum so that terms are real and positive. 

Analytical work in color superconductivity^ has demonstrated a rich 
phase structure at high density, and stimulated interest in QCD at non- 
zero baryon density. Several experiments have been proposed to probe 
matter at density of a few times nuclear matter density^. Even rudimen- 
tary information about the behavior of dense matter would be useful to the 
experimental program, as well as to the study of compact astrophysical ob- 
jects such as neutron stars. Recently^, we showed that QCD near a Fermi 
surface has positive, semi-definite measure. In the limit of low energies, 
the contribution of the remaining modes far from the Fermi surface can be 
systematically expanded, using a high density effective theory previously 
introduced by one of us^'^. This effective theory is sufficient to study phe- 
nomena like color superconductivity, although quantities like the equation 
of state may be largely determined by dynamics deep in the Fermi sea. 

The expansion about the Fermi surface is in powers of where fj, is 
the chemical potential. For this expansion to be controlled, the ultraviolet 
cutoff of our effective theory must be less than /i, or equivalently the scale 
of the physics of interest must be small relative to the chemical potential. 
In QCD at asymptotic density, the superconducting gap is exponentially 
small, so this condition is satisfied. However, it is also quite possible that 
at intermediate densities (e.g., those inside a neutron star) the gap is some- 
what smaller than /z, providing us with an additional small dimensionless 
parameter. Even if this is not the case, the power expansion of the effec- 
tive theory is qualitatively different from the usual perturbation in a^, and 
therefore worth exploring. 

2. Example: (1+1) Dimensions and Beyond 

We begin with an example that illustrates the basic ideas in a simple setting. 
Consider the Euclidean (1+1) action of non-relativistic fermions interacting 
with a gauge field A 




(1) 
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where e(p) is the energy as a function of momentum (e.g. e{p) w ^ H ). 

In (1) and below wc may consider it as a function of the operator {—idx+A). 
The dispersion relation in the presence of the chemical potential ep is: 
E{p) = e{p) — €F, and a low energy mode must have momentum close 
to ±pf, where e{±pF) = £f- The Fermi surface in (1+1) dimensions is 
reduced to the two points p = ±pf- Near these points we have 

E{p±pf)^±VfP , (2) 

where vf = dE/dp\pp is the Fermi velocity. 

The action (1) is not obviously positive. In fact, the operator in brackets 
[ • • • ] clearly has Hermitian as well as anti-Hermitian components, and 
hence complex eigenvalues. 

Let us assume that the gauge field has small amplitude and is slowly 
varying relative to the scale pF- Wc will extract the slowly varying compo- 
nent of the fermion field to construct a low energy effective theory involv- 
ing quasiparticles and gauge fields. This effective theory will have positive, 
semi-definite determinant. 

First, we extract the quasiparticle modes (we suppress the spin index 
in what follows) 

i:{x,T)=^Le+'P^^ + ^Pne-'P^^ , (3) 

where the functions ipL,R are slowly varying. To simplify the action, we use 
the identity 

^±ip^x E{-idx + A) e^^P^^ ii^{x) ^ ± VF{-id:c + A)ii^{x) , (4) 
to obtain'^ 



'eflf 



-I 



drdx ip\^{-dr + + idx - A)ijjL + 'tpni-dr + i(l> - id^ + A)'iljn 

(5)" 

We can write this in a more familiar form by introducing the Euclidean 
(1+1) gamma matrices 70.1,2 , which are Hermitian and can be taken as 
7j = CTj where a are the Pauli matrices. Using iI>l,r = 5(1±72)V' we obtain 



5eff = j drdx tp^^{d^ + iA^)^' = j drdx tppip . 



(6) 



''For simplicity wc set vp = 1. Alternatively it could be absorbed in the definition of 
the spatial 71 as in QCD below. 
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Since the gamma matrices are Hcrmitian, and the operator (9^ + iA^) is 
anti-Hermitian, the operator Ip in (6) has purely imaginary eigenvalues. 
However, because 72 anticommutes with ]/), the eigenvalues come in conju- 
gate pairs: given Ip (p = Xcf), we have 

P {72(f>) = -72-^ 4' = -I2M = -A(720n) • 

Hence the determinant det.^ = H''^*'^ is real and positive semi-definite. 

Thus, by considering only the low-energy modes near the Fermi points 
of the original model (1), wc obtain an effective theory with desirable pos- 
itivity properties. Note that it is necessary that the interactions (in this 
case, the background gauge field A) not couple strongly the low-energy 
modes to fast modes which arc far from the Fermi points. This is a reason- 
able approximation in many physical situations, where it is the interactions 
among quasiparticles that are of primary interest. In what follows, we will 
apply this basic idea to more complex models such as QCD. 

It is straightforward to go beyond (H-l) dimensions. Consider an elec- 
tron system, described by 

L = i,^[i^t-e{p)]i' + ^J^i'^i^, (7) 

where e{p) is the electron energy, a function of momentum p. It is interesting 
to note that the non-relativistic system already has a sign problem even at 
the zero density, /U = 0, though the free case does not suffer this, thanks to 
the separation of variables. In fact, it is quite unusual to have a system like 
vacuum QCD which has no sign problem. In Euclidean space the electron 
determinant is 

M = -dr- e{p) + IJL. (8) 

The first term in operator (8) is anti-Hermitian, while the rest are Hermi- 
tian. Since there is no constant matrix P in the spin space that satisfies 
= P M P~^, it has a sign problem in general. 
Let us decompose the fermion momentum as 

p=PF + L (9) 

Again, the Fermi momentum is defined to be a momentum at which the 
energy equals to the chemical potential at zero temperature: fj. = e{pF), 
and the Fermi velocity is defined as 
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If we are interested in low energies, / 
modes to get an effective operator, 

Meft = -dr -vf-T, (11) 

which has complex eigenvalues. However, when we include the —vp sector, 
we have Mept{vf)Meft{—vf) < (i.e., has real negative eigenvalues), 
assuming e{p) = e(— p). We again see that the sign problem in the electron 
system is alleviated in the low-energy effective theory. 

3. QCD 

Let us recall why the measure of dense QCD is complex in Euclidean space. 
We use the following analytic continuation of the Dirac Lagrangian to Eu- 
clidean space: 

xo -ix%, Xe ; 70 ^ il, 7, i'ys . (12) 

The Euclidean gamma matrices satisfy 

7^^ =7^ , {7^,7^} = 2^"^ (13) 

The Dirac-conjugated field, "ip = V'Wi is mapped into a field, still denoted 

as -0, which is independent of -0 and transforms as tp'^ under 50(4). Then, 
the grand canonical partition function for QCD is 

Z{^)= j dA^det{M)e-^^^'^\ (14) 

where <S'(^^) is the positive semi-definite gauge action, and the Dirac op- 
erator 

M = 7^7^^ + M7l, (15) 

where De = ds +iAE 'is the analytic continuation of the covariant deriva- 
tive. The Hermitian conjugate of the Dirac operator is 

M^ = -^%D%+„^% . (16) 

The first term in (15) is anti-Hermitian, while the second is Hermitian, 
hence the generally complex eigenvalues. When = 0, the eigenvalues are 
purely imaginary, but come in conjugate pairs (A, A*) ^, so the resulting 
determinant is real and positive semi-definite: 

detM = ]jA*A>0 . (17) 

''As before, note that 75 anti-commutes with M, so if Mcj) = A</>, then M'y^cj) = — 75M</> = 
-A75</>. 



<^Pf, we may integrate out the fast 
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In what follows we investigate the positivity properties of an effective 
theory describing only modes near the Fermi surface. To low energy modes, 
the curvature of the Fermi surface is not evident, and the positivity of the 
usual Dirac sea determinant is recovered. 

A system of degenerate quarks with a net baryon number asymmetry is 
described by the QCD Lagrangian density with a chemical potential /i, 

iqcD = Vii^ V - If^.F"""' + l^i'loi^, (18) 

where the covariant derivative = + iA^ and we neglect the quark 
mass for simplicity. 

At high density {n Aqcd), due to asymptotic freedom the energy 
spectrum of quarks near the Fermi surface is approximately given by a free 
Dirac eigenvalue equation, 

{a ■ p- fi)tjj± = E±il;±, (19) 

where a = 707 and ip± denote the energy cigcnfunctions with eigenvalues 
E± = ± \p\, respectively. At low energy E < ^, the states ^+ near the 
Fermi surface, \p\ ^ /z, are easily excited but which correspond to the 
states in the Dirac sea, are completely decoupled due to the presence of the 
energy gap /x provided by the Fermi sea. Therefore the appropriate degrees 
of freedom at low energy consist of gluons and ip+ only. 

Now, we wish to construct an effective theory describing the dynamics 
of ip+ by integrating out modes whose energy is greater than /x. Consider a 
quark near the Fermi surface, whose momentum is close to ij.vf- Without 
loss of generality, we may decompose the momentum of a quark into a Fermi 
momentum and a residual momentum as 

= P'V^ + l^, (20) 

where = (0, vp)- Since the quark energy is given as 

E=-f,+^m+fj,)^+ii, (21) 

the residual momentum should satisfy (Z|| +/i)^ + Zj^ < with Z|| = vpl-vp 

and l± = I — 

To describe the small excitations of the quark with Fermi momentum, 
IJ-vp, we decompose the quark fields as 

^Pix) = e'"^^^-^ [V+(t?F, x) + i>-{vF, x)] , (22) 

where 

i;±{vF,x)=P±{vF)e-'^''^-^ij{x) with P±(^J^) = l^^l^. (23) 
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The quark Lagrangian in Eq. (18) then becomes 

+ ■4)- {vf , x)'y° (2/i + i£>|| ) V- {vf , x) 

+ ['iij-{vF,x)ipj_ilj+{vF,x)+h.c.] (24) 

where ^1^ = (7°, vfVf ■ 7), 11=7"- 7^*, 5|| = Vi'D^ with = (1, vf), 

At low energy, we integrate out all the "fast" modes i/'- ^ind derive 
the low energy effective Lagrangian by matching all the one-liglit-particle 
irreducible amplitudes containing gluons and tp+ in loop expansion. The 
effects of fast modes will appear in the quantum corrections to the couplings 
of low energy interactions. At tree-level, the matching is equivalent to 
eliminating i/j- in terms of equations of motion: 

(25) 

Therefore, the tree-level Lagrangian for ?/'+ becomes 

LOff = - ^V;+7°(^i)'V'+ + • • • , (26) 

where the ellipsis denotes terms with higher derivatives. 

Consider the first term in our effective Lagrangian, which when contin- 
ued to Euclidean space yields the operator 

Meft = 7f • D{A). (27) 

Meft is anti-Hermitian and it anti-commutes with 75 , so it leads to a positive 
semi-definite determinant. However, note that the Dirac operator is not well 
defined in the space of 'tjj+{vF,x) (for fixed vf), since it maps ip+{vF,x) 
into ip+{—VF,x): 

i% P+V = P-i% V- (28) 

Since P-{vf) = P+{—vf), iptlJ+{vF,x) are tjj+{—VF,x) modes, or fluctu- 
ations of a quark with momentum —fivp- 

Thus far we have considered the quark velocity as a parameter labelling 
different sectors of the quark field. This is similar to the approach of heavy 
quark effective theory (HQET)^, in which the velocity of the heavy charm 
or bottom quark is almost conserved due to the hierarchy of scales between 
the heavy quark mass and the QCD scale. However, this approach contains 
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an ambiguity often rcfcrrccl to as "reparameterization invariance" , related 
to the non-uniqueness of the decomposition (20) of quark momenta into a 
large and residual component. In the dense QCD case, two tpivp, x) modes 
whose values of vf are not very different may actually represent the same 
degrees of freedom of the original quark field. In what follows we give a 
different formulation which describes all velocity modes of the quark field, 
and is suitable for defining the quasiparticle determinant. 

First, a more precise definition of the breakup of the quark field into 
Fermi surface modes. Using the momentum operator in a position eigen- 
state basis: p = —id, we construct the Fermi velocity operator: 

which is Hermitian, and a unit vector. 

Using the velocity operator, we define the projection operators P± as 
before and break up the quark field as, ijj{x) = ip+{x) + ip^{x), with ip± = 
P±ip. By leaving v as an operator we can work in coordinate space without 
introducing the HQET-inspired velocity Fourier transform which introduces 
wf as a parameter. If we expand the quark field in the eigcnstatcs of 
the velocity operators, we recover the previous formalism with all Fermi 
velocities summed up. 

The leading low-energy part of the quark action is given by 

L+=^P_iv)it^~4 + tijo)P+{v)i' ■ (30) 

As before, we define the fields 4'+ to absorb the large Fermi momentum: 

V+(a;) = e-"''^-^P+{v)ij{x). (31) 

Let us denote the eigenvalue v obtained by acting on the field ip (which has 
momentum of order /x) as vi (or v "large"), whereas eigenvalues obtained 
by acting on the effective field theory modes ijj^ arc denoted (or v "resid- 
ual"). If the original quark mode had momentum p with \p\ > ji (i.e. was 
a particle), then vi and Vr are parallel, whereas if \p\ < fj, (as for a hole) 
then Vr and vi arc anti-parallel. In the first case, we have P+(w;) = P-^{vr) 
whereas in the second case P+{vi) = P-{vr). Thus, the residual modes tp+ 
can satisfy either of P±{vr)ip+ = -0+, depending on whether the original ip 
mode from which it was derived was a particle or a hole. In fact, -0+ modes 
can also satisfy either of P±{vi)^^ — since they can originate from ^jj 
modes with momentum ~ +fiv as well as —fiv (both are present in the 



(29) 
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original measure: DipD-ip). So, the funetional measure for -0+ modes con- 
tains all possible spinor functions - the only restriction is on the momenta: 
1^0 1, 1^1 < A, where A is the cutoff. 

In light of the ambiguity between vi and Vr, the equation tjj = 6+''*^'"'^+ 
must be modified to 

ip = exp {+iij.x ■ V a ■ v)il)+ = exp {+iij,x ■ Vr a - Vr) tlJ+ , (32) 

where the factor of a. ■ v.,. corrects the sign in the momentum shift if Vr 
and vi are anti-parallel. In general, any expression with two powers of v 
is unaffected by this ambiguity. For notational simplicity we define a local 
operator 

X = f,x-va-v = (33) 
Taking this into account, we obtain the following action: 

L+ =Vi+e-^^(i^-4 + /Z7o)e+^^V+ • (34) 

We treat the ^ term separately from ii^ + /i7o since the former docs not 
commute with X, while the latter does. Continuing to Euclidean space, and 
using the identity P_7^P+ = 7]!^+, we obtain 

L+=^+7|'(a^ + zA^^)V+ , (35) 

where 

A^l = e-'^ A" e+'^ , (36) 

and all 7 matrices are Euclidean. The term containing A cannot be fully 
simplified because [v,A] ^ 0. Physically, this is because the gauge field 

carries momentum and can deflect the quark velocity. The redefined ip-^- 
modes are functions only of the residual momenta 1, and the exponential 
factors in the A term reflect the fact that the gluon originally couples to 
the quark field ip, not the residual mode ip+. 
The kinetic term in (35) can be simplified to 

71* = (37) 

since v ■ dv ■ "/ = 9-7 . The action (35) is the most general dimension 
4 term with the rotational, gauge invariance'' and projection properties 



°If we simultaneously gauge transform and V'+ in (35) the result is invariant. There 
is a simple relation between the gauge transform of the + fields and that of the original 

fields: U-f-{x) = U{x)e''-^ . Of course, the momentum-space support of the + gauge 
transform must be limited to modes less than the cutoff A. 
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appropriate to quark quasiparticles. Therefore, it is a general consequence 
of any Fermi liquid description of quark-like excitations. 

The operator in (35) is anti-Hermitian and leads to a positive, semi- 
definite determinant since it anti-commutes with 75. The corrections given 
in (26) arc all Hcrmitian, so higher orders in the expansion may re- 
introduce complexity. The structure of the leading term plus corrections 
is anti-Hermitian plus Hermitian, just as in the original QCD Dirac La- 
grangian with chemical potential. 

By integrating out the fast modes, the Euclidean QCD partition func- 
tion can be rewritten as 

Z(/x) = j dA+ detMeff(A+)e-^''«(^+). (38) 

The leading terms in the effective action for gluons (these terms are gen- 
erated when wo match our effective theory, with energy cutoff A, to QCD) 
also contribute only real, positive terms to the partition function: 



where A± = A — Ay and the Debye screening mass is M = ,^Nf /{2TT^)gs^ . 
Note that Landau damping is due to softer quark modes which have not 
been integrated out, and therefore do not contribute to matching. 

Although the HDET only describes low-energy modes, it still contains 
Cooper pairing interactions. This is because Cooper pairing, in which the 
quasiparticles have nearly equal and opposite momenta, is induced by glu- 
onic interactions with small energy and momentum transfer. That is, al- 
though a gluon exchange (or other interaction) which causes a large angular 
deflection of a quasiparticle 

\P) - \P'} 

must involve a large momentum transfer, and hence is not part of the 
effective theory, a Cooper pairing interaction 

\p, -p) W, -p') 

only involves a small energy and momentum transfer, even if the angle 
between p and p' is large. Hence, it is described by the leading order 
interaction between soft gluons and quarks in the effective theory (35). 
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4. Lattice Simulation 

The goal of this section is to give a method for simulating QCD at finite 
density. We will consider a chemical potential n much larger than Aqcd 
throughout, and divide the functional integral over quark excitations into 
two parts: (I) modes within a shell of width A of the Fermi surface, and (II) 
modes which are further than A from the Fermi surface. We will assume 
the hierarchy 

/X >> A >> Aqcd • (40) 

The quark determinant in region (I) is well approximated by the deter- 
minant of the leading operator in high density effective theory (HDET) as 
long as the first inequality in (40) is satisfied. As discussed in the previous 
section, it is positive and real. 

Here we will show that the contributions to the effective action for the 
gauge field from quark modes in region (II) are small and vanish as the A 
grows large compared to Aqcd • 

First consider the theory in Minkowski space. The Dirac operator is 

M = ip+ /i7o (41) 
and the Dirac equation can be written as 

ic^oV' = Htj} 

with 

H = ia-d-n (43) 

a Hermitian operator. The break up into regions (I) and (II) proceeds 

naturally in terms of energy eigenvalues of H (or Iq in the HDET notation). 
The low-lying modes in region (I) are particle states with spatial momenta 
satisfying \p\ « fj,. 

The analytic continuation of region (I) to Euclidean space leads to the 
HDET determinant considered previously. 

Modes in region (II) all have large energy eigenvalues, at least as large as 
A. In considering their effect on physics at the scale Aqcd, we can integrate 
them out in favor of local operators suppressed by powers of Aqcd/A. 

To make this concrete, consider the effective action for gauge fields with 
field strengths i^^^ of order Aqcd • The quark contribution to this effective 
action is simply the logarithm of the determinant we wish to compute. It 
can be expanded diagrammatically in graphs with external gauge field lines 
connected to a single quark loop. Restricting to region (II), we require that 



(42) 
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the quark modes in the loop have large H eigenvalues. Evaluating such 
graphs leads only to operators which are local in the external fields Afj_{x). 

The resulting renormalizable (dimension 4) operator is the finite density 
equivalent of F^^, except that due to the breaking of Lorentz invariance 
it contains separate time- and space-like components which represent the 
contribution of high-energy modes to the renormalization of the coupling 
constant, and Debye screening. These eff'ects do not introduce a complex 
component when continued to Euclidean space. 

Higher dimension operators, which involve additional powers of F^i, or 
covariant derivatives are suppressed by the scale A. These are pre- 
sumably the source of complex terms introduced to the effective action. 
However, due to the 1/A suppression they are dominated by the contribu- 
tion from the low- lying modes in region (I), which is necessarily non-local, 
but real. 

The logarithm of the Euclidean quark determinant will have the form: 

IndetM ~ 0(//) + (non - local, real) + 0(-^) (local, complex) , (44) 

where the first term is the (real, constant) free energy of non-interacting 
quarks, the second term is from the positive determinant in region (I) and 
the last term is the suppresscid, complex contribution from region (H). Only 
the last two terms depend on the gauge field Afj_{x) 

On the lattice, one can use the dominant dependence of detM on the 
first and second terms to do importance sampling. In order to keep the com- 
plex higher dimension operators (last term in (44)) small, it is important 
that the gauge field strengths are kept smaller than A^. One can impose 
this condition by using two different lattice spacings, Ug for the gluons and 
adet for the quarks, with Og > adet- The determinant is calculated on the 
finer add lattice, and is a function of plaquettes which are obtained by 
interpolation from the plaquettes on the coarser Og lattice. Interpolation 
can be defined in a natural way, since each lattice link variable is an 
clement of the gauge group, and one can connect any two points 51,(72 on 
the group manifold in a linear fashion: g{t) = gi + t{g2 — gi) , < t < 1. 

5. Inequalities and Anomaly Matching 

Positivity of the measure allows for rigorous QCD inequalities at asymptotic 
density. For example, inequalities among masses of bound states can be 
obtained using bounds on bare quasiparticle propagators. One subtlety 
that arises is that a quark mass term does not lead to a quasiparticle gap 
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(the mass term just shifts the Fermi surface). Hence, for technical reasons 
the proof of non-breaking of vector symmetries'*^ must be modified. (Naive 
apphcation of the Vafa-Witten theorem would preclude the breaking of 
baryon number that is observed in the color-flavor-locked (CFL) phase^'). 
A quasiparticlc gap can be inserted by hand to regulate the bare propagator, 
but it will explicitly violate baryon number. However, following the logic of 
the Vafa-Witten proof, any symmetries which are preserved by the regulator 
gap cannot be broken spontaneously. One can, for example, still conclude 
that isospin symmetry is never spontaneously broken (although sec below 
for a related subtlety). In the case of three flavors, one can introduce a 
regulator d with the color and flavor structure of the CFL gap to show 
rigorously that none of the symmetries of the CFL phase are broken at 
asymptotic density. On the other hand, by applying anomaly matching 
conditions'^, we can prove that the SU{i)A symmetries are broken. We 
therefore conclude that the CFL phase is the true ground state for three 
light flavors at asymptotic density, a result that was first established by 
explicit calculation'^'®''^. 

To examine the long-distance behavior of the vector current, we note 
that its correlator in a given background gauge field A can be written as 

{j;{vF,x)Jl{vF,y))'' = -Tr'y^T-S^{x,y;d)'y,T'S^{y,x;d), 

where the SU{Nf) flavor current J^{vf,x) = ■tp+{vp,x)'-ff^T""tjj+{vF,x). 
The propagator with S'J7(3)v-invariant IR regulator d is given as 

1 f°° 

S^{x, y; d) = {x\ - \y) = | dr {x\ e--(-^^) \y) 
where with D = d + iA 

'D-V d 



M = 70 



rft D-V 



(45) 



Since the eigenvalues of M are bounded from below by d, we have 
1 f°° e~'^^ 

{Aj^\y) < j^dTe-''^^/W)^/W) = —^^/W)^My), (46) 

where R= \x — y\. The current correlators fall off rapidly as — > oo; 
d^+detMeff(A) e-^=« {j^{vF,x)J^{vF,y))^* 

/. „-2dR f 

{j^{vF,x)J^{vF,y)) ^ <^J^ \{x\x)\\{y\y)\i47) 



< 
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where we used the Schwartz inequaUty in the first inequality, since the mea- 
sure of the effective theory is now positive, and equation (46) in the second 
inequahty. The IR regulated vector currents do not create massless modes 
out of the vacuum or Fermi sea, which implies that there is no Nambu- 
Goldstonc mode in the SU{3)v channel. Therefore, for three massless 
flavors SU{3)v has to be unbroken as in CFL. The rigorous result provides 
a non-trivial check on explicit calculations, and applies to any system in 
which the quasiparticle dynamics have positive measure. The case with 
non-zero quark masses is complicated, and requires careful consideration of 
the order of limits^. 

6. Conclusion 

The low-energy physics of dense fermionic matter, ranging from quark mat- 
ter to electronic systems, is controlled by modes near the Fermi surface. An 
effective Lagrangian describing the low-energy modes can be given in a sys- 
tematic expansion in powers of the energy scale over the chemical potential. 
The leading term in this expansion has a simple form, and we have shown 
that it leads to a real, positive Euclidean path integral measure. 

This observation opens the door to importance sampling in Monte Carlo 
simulations of dense matter systems. The key requirement is that the in- 
teractions do not strongly couple the low-energy modes to modes far from 
the Fermi surface. QCD at high density satisfies this requirement, as do 
all asymptotically free models. Electronic systems in which the important 
interactions involve momentum transfer less than the Fermi energy are in 
this category, although some idealized models such as the Hubbard model 
are not. We have given some proposals for how the positive effective theory 
might be simulated numerically. Ultimately, we hope that actual practi- 
tioners will develop even more practical methods. 

Finally, positivity has analytical applications as well, since it allows the 
use of rigorous inequalities. In QCD we obtain restrictions on symmetry 
breaking patterns at high density. Similar restrictions can probably be 
obtained for electronic systems with suitable interactions. 
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